Phase unwrapping with a rapid opensource minimum spanning tree algorithm (ROMEO)

PURPOSE: To develop a rapid and accurate MRI phase-unwrapping technique for challenging phase topographies encountered at high magnetic fields, around metal implants, or postoperative cavities, which is sufficiently fast to be applied to large-group studies including Quantitative Susceptibility Mapping and functional MRI (with phase-based distortion correction). METHODS: The proposed path-following phase-unwrapping algorithm, ROMEO, estimates the coherence of the signal both in space-using MRI magnitude and phase information-and over time, assuming approximately linear temporal phase evolution. This information is combined to form a quality map that guides the unwrapping along a 3D path through the object using a computationally efficient minimum spanning tree algorithm. ROMEO was tested against the two most commonly used exact phase-unwrapping methods, PRELUDE and BEST PATH, in simulated topographies and at several field strengths: in 3T and 7T in vivo human head images and 9.4T ex vivo rat head images. RESULTS: ROMEO was more reliable than PRELUDE and BEST PATH, yielding unwrapping results with excellent temporal stability for multi-echo or multi-time-point data. It does not require image masking and delivers results within seconds, even in large, highly wrapped multi-echo data sets (eg, 9 seconds for a 7T head data set with 31 echoes and a 208 × 208 × 96 matrix size). CONCLUSION: Overall, ROMEO was both faster and more accurate than PRELUDE and BEST PATH, delivering exact results within seconds, which is well below typical image acquisition times, enabling potential on-console application

Mitigating susceptibility-induced distortions in high-resolution 3DEPI fMRI at 7T

Geometric distortion is a major limiting factor for spatial specificity in high-resolution fMRI using EPI readouts and is exacerbated at higher field strengths due to increased B0 field inhomogeneity. Prominent correction schemes are based on B0 field-mapping or acquiring reverse phase-encoded (reversed-PE) data. However, to date, comparisons of these techniques in the context of fMRI have only been performed on 2DEPI data, either at lower field or lower resolution. In this study, we investigate distortion compensation in the context of sub-millimetre 3DEPI data at 7T. B0 field-mapping and reversed-PE distortion correction techniques were applied to both partial coverage BOLD-weighted and whole brain MT-weighted 3DEPI data with matched distortion. Qualitative assessment showed overall improvement in cortical alignment for both correction techniques in both 3DEPI fMRI and whole-brain MT-3DEPI datasets. The distortion-corrected MT-3DEPI images were quantitatively evaluated by comparing cortical alignment with an anatomical reference using dice coefficient (DC) and correlation ratio (CR) measures. These showed that B0 field-mapping and reversed-PE methods both improved correspondence between the MT-3DEPI and anatomical data, with more substantial improvements consistently obtained using the reversed-PE approach. Regional analyses demonstrated that the largest benefit of distortion correction, and in particular of the reversed-PE approach, occurred in frontal and temporal regions where susceptibility-induced distortions are known to be greatest, but had not led to complete signal dropout. In conclusion, distortion correction based on reversed-PE data has shown the greater capacity for achieving faithful alignment with anatomical data in the context of high-resolution fMRI at 7T using 3DEPI

Real‐time motion and retrospective coil sensitivity correction for CEST using volumetric navigators (vNavs) at 7T

PURPOSE: To explore the impact of temporal motion-induced coil sensitivity changes on CEST-MRI at 7T and its correction using interleaved volumetric EPI navigators, which are applied for real-time motion correction. METHODS: Five healthy volunteers were scanned via CEST. A 4-fold correction pipeline allowed the mitigation of (1) motion, (2) motion-induced coil sensitivity variations, Δ B 1 - , (3) motion-induced static magnetic field inhomogeneities, ΔB0 , and (4) spatially varying transmit RF field fluctuations, ΔB 1 + . Four CEST measurements were performed per session. For the first 2, motion correction was turned OFF and then ON in absence of voluntary motion, whereas in the other 2 controlled head rotations were performed. During post-processing Δ B 1 - was removed additionally for the motion-corrected cases, resulting in a total of 6 scenarios to be compared. In all cases, retrospective ∆B0 and - ΔB 1 + corrections were performed to compute artifact-free magnetization transfer ratio maps with asymmetric analysis (MTRasym ). RESULTS: Dynamic Δ B 1 - correction successfully mitigated signal deviations caused by head motion. In 2 frontal lobe regions of volunteer 4, induced relative signal errors of 10.9% and 3.9% were reduced to 1.1% and 1.0% after correction. In the right frontal lobe, the motion-corrected MTRasym contrast deviated 0.92%, 1.21%, and 2.97% relative to the static case for Δω = 1, 2, 3 ± 0.25 ppm. The additional application of Δ B 1 - correction reduced these deviations to 0.10%, 0.14%, and 0.42%. The fully corrected MTRasym values were highly consistent between measurements with and without intended head rotations. CONCLUSION: Temporal Δ B 1 - cause significant CEST quantification bias. The presented correction pipeline including the proposed retrospective Δ B 1 - correction significantly reduced motion-related artifacts on CEST-MRI

A method for the dynamic correction of B0-related distortions in single-echo EPI at 7 T

We propose a method to calculate field maps from the phase of each EPI in an fMRI time series. These field maps can be used to correct the corresponding magnitude images for distortion caused by inhomogeneity in the static magnetic field. In contrast to conventional static distortion correction, in which one 'snapshot’ field map is applied to all subsequent fMRI time points, our method also captures dynamic changes to B0which arise due to motion and respiration. The approach is based on the assumption that the non-B0-related contribution to the phase measured by each radio-frequency coil, which is dominated by the coil sensitivity, is stable over time and can therefore be removed to yield a field map from EPI. Our solution addresses imaging with multi-channel coils at ultra-high field (7 T), where phase offsets vary rapidly in space, phase processing is non-trivial and distortions are comparatively large. We propose using dual-echo gradient echo reference scan for the phase offset calculation, which yields estimates with high signal-to-noise ratio. An extrapolation method is proposed which yields reliable estimates for phase offsets even where motion is large and a tailored phase unwrapping procedure for EPI is suggested which gives robust results in regions with disconnected tissue or strong signal decay. Phase offsets are shown to be stable during long measurements (40 min) and for large head motions. The dynamic distortion correction proposed here is found to work accurately in the presence of large motion (up to 8.1°), whereas a conventional method based on single field map fails to correct or even introduces distortions (up to 11.2 mm). Finally, we show that dynamic unwarping increases the temporal stability of EPI in the presence of motion. Our approach can be applied to any EPI measurements without the need for sequence modification

In vivo phase imaging of human epiphyseal cartilage at 7 T

PURPOSE: To assess the potential clinical utility of in vivo susceptibility-weighted imaging and quantitative susceptibility mapping of growth cartilage in the juvenile human knee at 7 T. METHODS: High-resolution gradient-echo images of the knees of six healthy children and adolescents aged 6 to 15 were acquired with a 28-channel coil at 7 T. Phase images from the coils were combined using a short echo-time reference scan method (COMPOSER). RESULTS: Veins oriented perpendicular to the static B0field appeared doubled in susceptibility-weighted imaging, but not quantitative susceptibility mapping. Veins and layers in the cartilage were visible in all children up to the age of 13. CONCLUSION: Phase imaging using susceptibility-weighted imaging and quantitative susceptibility mapping allows the in vivo visualization of veins and layers in human growth cartilage

Key clinical benefits of neuroimaging at 7 T

The growing interest in ultra-high field MRI, with more than 35.000 MR examinations already performed at 7 T, is related to improved clinical results with regard to morphological as well as functional and metabolic capabilities. Since the signal-to-noise ratio increases with the field strength of the MR scanner, the most evident application at 7 T is to gain higher spatial resolution in the brain compared to 3 T. Of specific clinical interest for neuro applications is the cerebral cortex at 7 T, for the detection of changes in cortical structure, like the visualization of cortical microinfarcts and cortical plaques in Multiple Sclerosis. In imaging of the hippocampus, even subfields of the internal hippocampal anatomy and pathology may be visualized with excellent spatial resolution. Using Susceptibility Weighted Imaging, the plaque-vessel relationship and iron accumulations in Multiple Sclerosis can be visualized, which may provide a prognostic factor of disease. Vascular imaging is a highly promising field for 7 T which is dealt with in a separate dedicated article in this special issue. The static and dynamic blood oxygenation level-dependent contrast also increases with the field strength, which significantly improves the accuracy of pre-surgical evaluation of vital brain areas before tumor removal. Improvement in acquisition and hardware technology have also resulted in an increasing number of MR spectroscopic imaging studies in patients at 7 T. More recent parallel imaging and short-TR acquisition approaches have overcome the limitations of scan time and spatial resolution, thereby allowing imaging matrix sizes of up to 128×128. The benefits of these acquisition approaches for investigation of brain tumors and Multiple Sclerosis have been shown recently. Together, these possibilities demonstrate the feasibility and advantages of conducting routine diagnostic imaging and clinical research at 7 T

The Impact of Echo Time Shifts and Temporal Signal Fluctuations on BOLD Sensitivity in Presurgical Planning at 7 T.

OBJECTIVES: Gradients in the static magnetic field caused by tissues with differing magnetic susceptibilities lead to regional variations in the effective echo time, which modifies both image signal and BOLD sensitivity. Local echo time changes are not considered in the most commonly used metric for BOLD sensitivity, temporal signal-to-noise ratio (tSNR), but may be significant, particularly at ultrahigh field close to air cavities (such as the sinuses and ear canals) and near gross brain pathologies and postoperative sites. MATERIALS AND METHODS: We have studied the effect of local variations in echo time and tSNR on BOLD sensitivity in 3 healthy volunteers and 11 patients with tumors, postoperative cavities, and venous malformations at 7 T. Temporal signal-to-noise ratio was estimated from a 5-minute run of resting state echo planar imaging with a nominal echo time of 22 milliseconds. Maps of local echo time were derived from the phase of a multiecho GE scan. One healthy volunteer performed 10 runs of a breath-hold task. The t-map from this experiment served as a criterion standard BOLD sensitivity measure. Two runs of a less demanding breath-hold paradigm were used for patients. RESULTS: In all subjects, a strong reduction in the echo time (from 22 milliseconds to around 11 milliseconds) was found close to the ear canals and sinuses. These regions were characterized by high tSNR but low t-values in breath-hold t-maps. In some patients, regions of particular interest in presurgical planning were affected by reductions in the echo time to approximately 13-15 milliseconds. These included the primary motor cortex, Broca's area, and auditory cortex. These regions were characterized by high tSNR values (70 and above). Breath-hold results were corrupted by strong motion artifacts in all patients. CONCLUSIONS: Criterion standard BOLD sensitivity estimation using hypercapnic experiments is challenging, especially in patient populations. Taking into consideration the tSNR, commonly used for BOLD sensitivity estimation, but ignoring local reductions in the echo time (eg, from 22 to 11 milliseconds), would erroneously suggest functional sensitivity sufficient to map BOLD signal changes. It is therefore important to consider both local variations in the echo time and temporal variations in signal, using the product metric of these two indices for instance. This should ensure a reliable estimation of BOLD sensitivity and to facilitate the identification of potential false-negative results. This is particularly true at high fields, such as 7 T and in patients with large pathologies and postoperative cavities

In vivo MRI of the human finger at 7 T

PURPOSE: To demonstrate a dedicated setup for ultrahigh resolution MR imaging of the human finger in vivo. METHODS: A radiofrequency coil was designed for optimized signal homogeneity and sensitivity in the finger at ultrahigh magnetic field strength (7 T), providing high measurement sensitivity. Imaging sequences (2D turbo-spin echo (TSE) and 3D magnetization-prepared rapid acquisition gradient echo (MPRAGE)) were adapted for high spatial resolution and good contrast of different tissues in the finger, while keeping acquisition time below 10 minutes. Data was postprocessed to display finger structures in three dimensions. RESULTS: 3D MPRAGE data with isotropic resolution of 200 µm, along with 2D TSE images with in-plane resolutions of 58 × 78 µm2 and 100 × 97 µm2 , allowed clear identification of various anatomical features such as bone and bone marrow, tendons and annular ligaments, cartilage, arteries and veins, nerves, and Pacinian corpuscles. CONCLUSION: Using this dedicated finger coil at 7 T, together with adapted acquisition sequences, it is possible to depict the internal structures of the human finger in vivo within patient-compatible measurement time. It may serve as a tool for diagnosis and treatment monitoring in pathologies ranging from inflammatory or erosive joint diseases to injuries of tendons and ligaments to nervous or vascular disorders in the finger. Magn Reson Med 79:588-592, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited

A comparison of static and dynamic ∆B0 mapping methods for correction of CEST MRI in the presence of temporal B0 field variations

Purpose: To assess the performance, in the presence of scanner instabilities, of three dynamic correction methods which integrate ∆B 0 mapping into the chemical exchange saturation transfer (CEST) measurement and three established static ∆B 0 -correction approaches. Methods: A homogeneous phantom and five healthy volunteers were scanned with a CEST sequence at 7 T. The in vivo measurements were performed twice: first with unaltered system frequency and again applying frequency shifts during the CEST acquisition. In all cases, retrospective voxel-wise ∆B 0 -correction was performed using one intrinsic and two extrinsic [prescans with dual-echo gradient-echo and water saturation shift referencing (WASSR)] static approaches. These were compared with two intrinsic [using phase data directly generated by single-echo or double-echo GRE (gradient-echo) CEST readout (CEST-GRE-2TE)] and one extrinsic [phase from interleaved dual-echo EPI (echo planar imaging) navigator (NAV-EPI-2TE)] dynamic ∆B 0 -correction approaches [allowing correction of each Z-spectral point before magnetization transfer ratio asymmetry (MTR asym) analysis]. Results: All three dynamic methods successfully mapped the induced drift. The intrinsic approaches were affected by the CEST labeling near water (∆ω < |0.3| ppm). The MTR asym contrast was distorted by the frequency drift in the brain by up to 0.21%/Hz when static ∆B 0 -corrections were applied, whereas the dynamic ∆B 0 corrections reduced this to <0.01%/Hz without the need of external scans. The CEST-GRE-2TE and NAV-EPI-2TE resulted in highly consistent MTR asym values with/without drift for all subjects. Conclusion: Reliable correction of scanner instabilities is essential to establish clinical CEST MRI. The three dynamic approaches presented improved the ∆B 0 -correction performance significantly in the presence of frequency drift compared to established static methods. Among them, the self-corrected CEST-GRE-2TE was the most accurate and straightforward to implement

The clinical relevance of distortion correction in presurgical fMRI at 7 T

Presurgical planning with fMRI benefits from increased reliability and the possibility to reduce measurement time introduced by using ultra-high field. Echo-planar imaging suffers, however, from geometric distortions which scale with field strength and potentially give rise to clinically significant displacement of functional activation. We evaluate the effectiveness of a dynamic distortion correction (DDC) method based on unmodified single-echo EPI in the context of simulated presurgical planning fMRI at 7 T and compare it with static distortion correction (SDC). The extent of distortion in EPI and activation shifts are investigated in a group of eleven patients with a range of neuropathologies who performed a motor task. The consequences of neglecting to correct images for susceptibility-induced distortions are assessed in a clinical context. It was possible to generate time series of EPI-based field maps which were free of artifacts in the eloquent brain areas relevant to presurgical fMRI, despite the presence of signal dropouts caused by pathologies and post-operative sites. Distortions of up to 5.1 mm were observed in the primary motor cortex in raw EPI. These were accurately corrected with DDC and slightly less accurately with SDC. The dynamic nature of distortions in UHF clinical fMRI was demonstrated via investigation of temporal variation in voxel shift maps, confirming the potential inadequacy of SDC based on a single reference field map, particularly in the vicinity of pathologies or in the presence of motion. In two patients, the distortion correction was potentially clinically significant in that it might have affected the localization or interpretation of activation and could thereby have influenced the treatment plan. Distortion correction is shown to be effective and clinically relevant in presurgical planning at 7 T